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We theoretically derive the polarization-resolved intensity distribution of a 
TM-polarized fundamental Gaussian beam reflected by an air-glass plane 
interface at Brewster incidence. The reflected beam has both a dominant 
(TM) and a cross-polarized {TE) component, carried by a TEMio and a 
TEMqi Hermite-Gaussian spatial mode, respectively. Remarkably, we find 
that the TE'-mode power scales quadratically with the angular spread of the 
incident beam and it is comparable to the TM-mode power. Experimental 
confirmations of the theoretical results are also presented. © 2009 Optical 
Society of America 
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When a beam of light impinges upon a plane interface separating two transparent media, 
it produces reflected and transmitted beams. In 1815 the Scottish physicist David Brew- 
ster discovered the total polarization of the reflected beam at the angle 9b since named after 
him [1] . From his observations he was also able to empirically determine the celebrated equa- 
tion, known as Brewster's law, tan#g = ni/ri2, where ri\ and ri2 are the respective refractive 
indices of the two media. Several articles have been published on theory and experiments 
about Brewster's law for beams with non-planar wave fronts. Fainman and Shamir [2] ad- 
dressed reflection of an (isotropic) spherical wave, by a Brewster angle polarizer; they found 
a cross-polarized component. More recently, Kohazi-Kis [3] has theoretically derived and 
experimentally confirmed cross-polarization effects occurring at Brewster incidence. Shortly 
afterwards, Li and Vernon [4] addressed the same problem, using a microwave Gaussian 
beam; they do not mention cross polarization. Consequences of cross-polarization coupling 
(XPC) at dielectric interfaces were also theoretically investigated by Nasalski and coworkers 
in a series of interesting papers [5,6]. 

However, these studies fail to compare the intrinsic XPC due to the natural angular spread 
(namely, the focusing) of the incident beam, and the non-intrinsic one caused by reflection 
of such a beam at a dielectric interface. The main aim of this Letter is to fill this gap. 

The structure of this Letter is as follows. We first solve the general problem of the reflection 
of a polarized fundamental Gaussian beam at the plane interface between two optical media 
[7,8]. Next, we derive analytical expressions for the polarization-dependent transverse spatial 
profiles of the reflected beam. From this result, we are able to prove that non-intrinsic 
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XPC scales quadratically with the angular spread 9 of the incident beam, as opposed to 
the intrinsic XPC that scales as the fourth power of 9q. Finally, we present experimental 
confirmations of our theoretical findings. 



Fig. 1. (Color online) Geometry of beam reflection at the air-medium interface. 9b is the 
Brewster angle. 

Consider a monochromatic beam of light incident upon a plane interface that separates air 
from glass. With n = n a i r /n g i ass we denote the ratio between the two refractive indices. As 
the beam meets the interface coming from the air side, it will be convenient to take the axis 
z of the laboratory Cartesian frame K = (O, x, y, z) normal to the interface and directed 
from the air to the glass. Moreover, we choose the origin O in a manner that the plane 
interface has equation z = 0. The air-glass interface, the incident and the reflected beams 
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are pictorially illustrated in Fig. 1. In addition to the laboratory frame, we use a Cartesian 
frame Ki = (O, Xi, jji, z{) attached to the incident beam and another one K r = (O, x r , y r , z r ) 
attached to the reflected beam. Let k = k^Zi and k denote the central and noncentral wave 
vectors of the incident beam, respectively, with |k| = |ko| = ko- Then, the electric field of the 
incident beam can be written as a linear superposition of the fundamental vector plane-wave 
mode functions Xx{k) = eA(k)exp(ik ■ r) with complex amplitudes a\(k), as follows: 

E'(r)=£ I a A (k)*A(k)d 2 fc T , (1) 

A=l 

where ky = k — k (k ■ k)//cQ is the transverse part of k, and we choose the polarization unit 
basis vectors as ei(k) = e 2 (k) x k/fc , and e 2 (k) = zxk/|zxk| [9]. Here z is a real unit 
vector directed along the laboratory axis z. By definition, ei(k) lies in the plane of incidence 
containing both the wave vector k and z, while e 2 (k) is orthogonal to such a plane. A plane 
wave whose electric field vector is parallel to either ei(k) or e 2 (k), is referred to as either 
a TM or a TE wave, respectively. The symbols S for TE and P for TM, are also widely 
used. In Eq. ([T|) a\(k) = A(k) a\(k), where A(k) and a\(k) are the scalar and the vector 
spectral amplitudes of the field, respectively. Here we consider a monochromatic Gaussian 
beam, whose spectral amplitude A(k) is localized in k space, centered at the central wave 
vector k = k z,i, on the sphere of equation cu 2 (k) = c 2 /cq, namely 

A{k) = e-^^e^ 1 -!^! 2 ) 172 , (2) 

where 8 = 2/(k wo) is the diffraction-defined angular aperture of the incident beam [10] 
which has, by hypothesis, a minimum diameter (spot size) equal to 2u> located at Zi = —d. 
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The vector spectral amplitudes are defined as a\(k) = e A (k) • /, where / = (/pXj + fsfi), 
with \fp\ 2 + I/5 1 2 = 1, is a complex- valued unit vector that fixes the polarization of the 
incident beam. 

When the latter is reflected at the interface, each vector mode function changes according 

to 

Xx(k)^r x (k)xx(k), (3) 

where r^k) and r 2 (k) are the Fresnel reflection amplitudes for TM and TE waves, respec- 
tively [11], and k = k — 2 z (z ■ k) is sets by the law of specular reflection [12]. If we substitute 
Eq. ([3]) into Eq. (TT|), we obtain 

2 r 

E'(r) 1 ► E*(r) = £ / a A (k)r A (k)x A (k) d 2 £; T , (4) 

A=l •* 

where k = koz r , by definition. The expression for the magnetic field B- R (r) of the re- 
flected beam may be obtained from the equation above via the straightforward substitutions 
a A (k) — > fe A (k)/c, where foi(k) = -a 2 (k)r 2 (k)/ri(k), and b 2 (k) = ai(k)r 1 (k)/r 2 (k). 

From Eq. (|2j) it follows that A(k) ~ for those wave vectors k lying outside the paraxial 
domain V = {k : |ky/fco| < ^0}, with 6 1 for well-collimated beams. This allows us to 
calculate analytically E R (and, similarly, B^) via a power series expansion for the integrand 
of Eq. (j3J) about the point k = ko up to and including second order terms in kr/^o- In 
practice, we extend to the problem at hand, the perturbative approach introduced by Lax 
et al, [13], and further developed by Deutsch and Garrison [14]. It is easy to see that if with 
7(k) = arccos(k • k //c 2 ) we denote the angle between the central wave vector k and the 
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non-central one k, then we can write |kr//c | = shi7(k) ~ 7(k), where 7(k) < 6q <C 1, being 
8q the natural small parameter for the power series expansion [14]. Explicit expressions for 
the power series expansions of both E i? (r) and B i? (r) are given in Appendix A. 

From the knowledge of both E R (r) and B^r) it is possible to calculate the intensity spatial 
distribution (i.e., the beam profile) I(r) of the reflected beam, as the flux of the cycle-averaged 
Poynting vector S oc Re(E fi x B^*) across a surface perpendicular to the central direction 
of propagation z r , namely I(r) oc S • z r . Since (E R x B fl *) • z r = E^B^* — Ey r B R r *, we 
can write /(r) = I P (r) + I s (r), where I P (r) = Re(E*B**), and I s (r) = Re(-E^B^*) are 
the intensity distributions produced by the P- and S'-polarized components of the reflected 
beam, respectively. Here P and S polarization directions are defined with respect to the 
central plane of incidence containing z, ko, and ko. After a lengthy but straightforward 
calculation it is not difficult to obtain, for a P-polarized incident beam (i.e., for the choice 
fp = l and fs = 0), the following power series expansions: 

Jp(r)// (r) =r\ + 6,uX + 6 2 (v + pX 2 + qY 2 ) , (5) 
I s (r)/I (r)=e 2 sY 2 , (6) 

where X = x r /w , Y = y r /w , Z = (z r + d)/L, and Jo( r ) = ex P {^^ X \+z^ /(-'■ + ^ s 
intensity distribution of the incident beam. In Eq. ([5]) the transverse coordinates x r and y r 
are normalized with respect to the beam waist wq, while the longitudinal coordinate z r is 
normalized with respect to the Raleigh range L = kw"l/2 of the beam. In Eqs. ([5]) and ([6]), r P 
and r$ are the Fresnel reflection amplitude for P and S waves, respectively, evaluated at the 
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central angle of incidence 9 = arccos(k • z/k ), while u, v, p, q and s are some complicated 
functions of Z, 9, rp,r$ and their derivatives [15] whose explicit form is given in Appendix 
B. 

However, at the Brewster angle of incidence 9 = 9p = arctan(n), only p and s take a 
non-zero value, namely 

p(l + Z 2 ) = {dr P /d9)\l , ,(1 + Z 2 ) = rl/n\ B , (7) 

and Eqs. (J5j) and ([6]) reduces to 

Jp(r)// (r) = 9\pX\ / 5 (r)// (r) = 9 2 sY 2 , (8) 

respectively. From this result it immediately follows that the ratio p between the power of 
the S and the P components of the reflected beam at Brewster incidence, is simply equal to 

s/p, 

ni s (r)AXdY fr a 1 V 
P {{I P (r)dXdY Kn&rp/dO),^ { > 

This simple result is remarkable since it shows that p is independent from the waist w of 
the incident beam, that is, p take the same value for either a well-collimated or a strongly- 
focussed beam. 

Equations ([MS]) represent the main theoretical result of this Letter. In particular, Eq. ([6]) 
shows that non-intrinsic XPC generates an intensity Is(r) that scales quadratically with 9 . 
This behavior cannot be ascribed to the intrinsic XPC exhibited by the incident beam since 
the latter is due to the beam-divergence only and the consequent cross-polarized intensity 
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scales with 9^ [16,17]. Moreover, Eq. (jSJ) shows two important things. First, we see that 
even at Brewster angle of incidence the extinction of a P-polarized beam is not perfect, 
as (drp/d9)\ g ^ and rg\ e ^ 0. Thus, although the input beam is P-polarized, a S 
component appears after reflection. Second, the two expression in Eq. (jSJ) show that after 
reflection the cylindrical symmetry about the axis of propagation of the beam is lost and 
two orthogonally polarized TEMio and TEMoi modes are generated. The fact that Is has a 
TEM i profile, as opposed to the cloverleaf TEMn pattern typical of the cross polarization 
intensity of the incident beam, is consistent with the hypothesis that such a term originates 
from non-intrinsic XPC and it is not a simple beam-divergence effect. 

We verified these theoretical results in our laboratory by using a Super-luminescent Light 
Emitting Diode (SLED) operating at A = 820 nm (InPhenix IPSDD0802) as a light source. 
The output of the SLED was first spatially filtered by a single-mode optical fiber to prepare 
the input beam into the fundamental Gaussian mode, and then collimated by a microscope 
objective to produce a very large beam waist, (wq = 1.64 mm) before passing across a 
polarizer selecting P-polarization. A second lens put behind the polarizer generated the 
desired waist for the input beam. The so-prepared beam was sent upon the surface of a right- 
angle BK7 glass (n = 1.51) prism mounted on a precision rotation stage with a resolution of 
9 x 10 _6 rad (Newport URS-BCC) to accurately determine the Brewster angle (9b — 56.49°). 
Finally, the polarization-dependent beam intensity profiles after reflection were recorded 
by a CCD-based Beam Intensity Profiler (Spiricon LBA-FW-SCOR-20) mounted behind a 
polarizer put along the axis of the reflected beam, at large distance from the interface (far 
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field measurement). 



Theory Experiment 




Fig. 2. (Color online) Calculated and measured intensity transverse spatial profiles of the 
P- and S'-polarized modes of the reflected beam. The beam waist of the incident beam was 
wo = 34 (Ltm. 

A qualitative comparison between calculated and measured intensity distributions is shown in 
Fig. 2. The measured ratio p exp of the ^-polarization component power to the P-polarization 
one at Brewster incidence, was p exp = 0.20 ±0.05 in excellent agreement with the theoretical 
prediction of Eq. ([9]) giving, for BK7 glass, p t h = 4n 4 /(l + n 2 ) 4 = 0.18. 



9 



In conclusion, we found that when a TM-polarized fundamental Gaussian beam is re- 
flected at Brewster incidence it generates a two-mode beam with both a dominant and a 
cross-polarized component. The intensity of the latter scales quadratically with the angular 
divergence of the incident beam and can be, therefore, orders of magnitude bigger than the 
intrinsic cross-polarized intensity of the incident beam that scales with the fourth power of 



This project is supported by FOM. 
Appendix A 

As a result of the power series expansion truncated at second order terms in Q , we found the 
following expressions for the electric and magnetic fields for the reflected beam, evaluated in 
the beam frame K r : 



00- 




(10) 
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where X = x r /w , Y = y r /w Q , Z = (z r + d)/L, and 

rj, fl 2 j , (r P -2r P ) ,r P cotfl *(r P + r 5 ) cot 2 
-^x r = r P + #0^ " = + — < I— — — h «: 



z - Z 2 2(i-Z) 2(i - Z) 



i - Z 



X 



II 

2 VP 



r'p — 2r P ) 



(i - Z)= 



rpcot^ 2(r P + r 5 ) cot 2 6» 



L(i-z)= 



E yT =e ^ rp+ : s l coi0 + elxY 



E Zr =9 - 



i - Z 

rpX ■ el 



r'p + r' s ) cot — (r P + rg) esc 2 9 



i - Z 



+ i 



rp + rs) cot 9 



(i - zy 



(12) 
(13) 



2(i - Z) 2(i - Z) 



X* 



(i - zy 



y2 i r P + r s) cot 9 



(i - zy 



and 



(14) 



B I , = -e Y (rF + r f ot9 + elxY 

% — Z 



(r'p + r' s ) cot 9 (r P + r s ) cot 2 9 



B Vr =r P + 9 X 



(i - zy (i - zy 

r'p 91 j. r'p cot 9 .r" P + 2r s - 2(r P + r s ) esc 2 9 
i-Z + 1 2\2(i-Z) +% 2ji - Z) 



+ X 



(r'p - r P ) 

(i - zy 



Y 



r'p cot 9 rp + 2r$ — 2(rp + rs) esc 2 9 



zy 



zy 



9 n Y-^ + 9 2 n XY 



B Zr -U Q 1 T u {y 

% — z 



zy 



(r P + rs) cot 9 

(i - zy 



(15) 



(16) 
(17) 



The expressions above have been calculated for a P-polarized incident beam, and rp and rs 
are the Fresnel reflection amplitude for P and S waves, respectively, evaluated at the central 
angle of incidence 9 = arccos(k ■ &/k ). Moreover, we have used the notation r' A = dr A /d9 
and r" A = d 2 rA/d9 2 , where A £ {P, S}. It is worth noting that the cross-polarization terms 
at first order in 9q present in the expressions of E Vr and B Xr , have the same functional 
dependence ~ (rp + rs) given in Eqs. (16) and (19) of Ref. [6]. Therefore, since for the 
geometrically reflected beam (which is the specular image of the incident one) we have 
r P = 1, r s = —1, it is clear that first order cross-polarization terms disappear. 
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Appendix B 

In Eqs. ([MS]) u, v, p, q and s are some functions of Z, 9, r P ,r s and their derivatives. Their 
explicit form is given below: 

Z 



r P [r' P + 2r s + r' P cot 9 - 2 (r p + r s ) esc 2 9] 

(19) 



2(1 + Z 2 ) 

_ -r P (l-Z 2 ) (2^-3r P ) + 2r^(l + Z 2 ) 
P- 4(1 + Z 2 ) 2 ' [Z[)) 

_ r P (l - ^ 2 ) esc 2 [5r P + 4r 5 + (3r P + 4r s ) cos 29 - 2r' P sin 20] 
q ~ 2(1 + Z 2 ) 2 ' ( } 

_ (rp + r s ) 2 cot 2 

s _ i + z 2 • {ZZ) 

These functions look quite complicated, however u, v and q are proportional to rp so that 
they disappear when evaluated at Brewster angle of incidence. In passing, we mention here 
that the linear term 9quX in Eq. (jHJ) generates the Goos-Hanchen shift in the reflected beam. 
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